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We present the first measurement of elliptic ( 02 ) and triangular ( 113 ) flow in high-multiplicity 
3 He+Au collisions at ^/s NN = 200 GeV. Two-particle correlations, where the particles have a large 
separation in pseudorapidity, are compared in 3 He+Au and in p+p collisions and indicate that 
collective effects dominate the second and third Fourier components for the correlations observed in 
the 3 He+Au system. The collective behavior is quantified in terms of elliptic V 2 and triangular V 3 
anisotropy coefficients measured with respect to their corresponding event planes. The V 2 values are 
comparable to those previously measured in d+Au collisions at the same nucleon-nucleon center- 
of-mass energy. Comparison with various theoretical predictions are made, including to models 
where the hot spots created by the impact of the three 3 He nucleons on the Au nucleus expand 
hydrodynamically to generate the triangular flow. The agreement of these models with data may 
indicate the formation of low-viscosity quark-gluon plasma even in these small collision systems. 

PACS numbers: 25.75.Dw 


The study of high-energy heavy ion collisions at the 
Relativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider has produced abundant evidence for the 
formation of Quark-Gluon Plasma (QGP). Observation 
of strong elliptic and triangular flow in these A+A col¬ 
lisions indicates that the QGP has very small viscosity 
and behaves like a nearly perfect fluid m- Recently, 
measurements at the large hadron collider in very high 
multiplicity events from collisions of p+p and p+Pb have 
revealed similar particle emission patterns gHS] . Such 
features have also been detected in d+Au collisions at 
RHIC [QUIT]. Explanations of the data in terms of the 
formation of small droplets of QGP which then expand 
hydrodynamically m compete with alternatives involv¬ 
ing novel initial-state effects (e.g. glasma models H31). 
RHIC is uniquely suited to test these competing theories 
by its ability to engineer the size and shape of the ini¬ 
tial reaction zone through collisions of p+Au, d+Au, and 
3 He+Au, as proposed in Ref. m ■ In the case when small 
QGP droplets are formed, the latter two systems should 
have the strongest elliptic and triangular flow patterns, 
respectively. 

The azimuthal anisotropy of produced particles can 
be quantified by the Fourier coefficients v n in the ex¬ 
pansion of the particles’ distribution as: dN/d(j) oc 1 + 
^ n=1 2 v n cos(n(</>— 4'„)) |T5] , where n is the order of the 
harmonic, <fi is the azimuthal angle of particles of a given 
type, and \I/ n is the azimuthal angle of the nth-order 
event plane. In this Letter, the elliptic ( V 2 ) and triangu¬ 
lar ( 1 ) 3 ) flow for inclusive charged hadrons produced at 
midrapidity |?y| < 0.35 in high-multiplicity 3 He+Au colli¬ 
sions at ,ys NN = 200 GeV are measured in the PHENIX 
experiment with respect to ’L 2 and ^3 event planes. The 
signatures of collective motion beyond known nonflow 
correlation effects (e.g. jets, resonances, etc.) are also 
examined by comparing long-range correlations between 
3 He+Au and p+p collisions. 

A full description of the PHENIX experimental setup 
is given in Ref. [16]. Charged particles are reconstructed 
in the two PHENIX central-arm tracking systems com¬ 


prising drift chambers and multi-wire proportional pad 
chambers (PC) 17], Each arm covers ir/2 in azimuth 
and \p\ < 0.35. The drift-chamber tracks are matched to 
hits in the outermost PC layer (PC3), reducing the con¬ 
tribution of tracks originating from decays and photon 
conversions. 

The observed event-plane angles v k® bs are measured in 
different pseudorapidity ranges by the beam-beam coun¬ 
ters [THJ (BBC) and forward-silicon-vertex [15] (FVTX) 
detectors. The PHENIX experiment has two BBCs, each 
comprising 64 quartz Cerenkov radiators read out by 
photomultiplier tubes (PMTs), subtending pseudorapid¬ 
ity (3.0 < |ry| < 3.9). The FVTX [19] detector comprises 
two identical endcap assemblies, located symmetrically 
in the north and south directions. Charged particles can 
be detected with a high efficiency(> 95%) using a clus¬ 
ter of mini-strip hits. The event planes are measured 
by the BBC in the Au-going (south) direction (BBC-S), 
which covers —3.9 < p < —3.0, and by the reconstructed 
clusters in the FVTX in the Au-going (south) endcap 
(FVTX-S), which covers —3.0 < p < —1.0. 

The 3 He+Au data for this analysis were obtained in 
the 2014 run of the PHENIX experiment and include 
1.6 billion minimum-bias (MB) triggered events and 480 
million high-multiplicity (HM) triggered events. The MB 
trigger is defined as a coincidence between the north and 
south BBCs requiring one or more photomultiplier tubes 
firing in each, capturing 88±4% of the total inelastic 
3 He+Au cross section. The HM trigger is based on the 
MB trigger, but additionally requires more than 48 pho¬ 
tomultiplier tubes firing in the BBC-S. 

The event centrality class in 3 He+Au collisions is de¬ 
termined as a percentile of the total charge (%)Q BBC ~ S 
measured in the BBC-S [TSJ [50H22I- The distribution 
of BBC-S PMT charge sum is shown in Fig. [T| for both 
MB and HM triggered events. The distribution for HM 
events has been scaled down by the relative online trigger 
prescale factor. For the MB sample, a threshold on the 
BBC-S charge sum shown in Fig. [l] is applied to select 
the top 5% central 3 He+Au collisions. In total, 400 mil- 
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lion events have been used for the measurement of V 2 and 
Vs including both HM and MB triggered events having 
BBC-S charge sum above the threshold. 

Events with high BBC-S sum charge have been sim¬ 
ulated with a Monte Carlo Glauber model [23J [Mj , fol¬ 
lowing the procedure detailed in Ref. |T5]. The average 
number of binary collisions (A c 0 n), participants (A part ) 
and initial-state eccentricities (£ 2 , £ 3 ) were found to be 
consistent, within uncertainties, between the MB with 
the 0%-5% central selection, and the HM events being 
used here. Note that in this Glauber calculation the spa¬ 
tial distribution from each participant is smeared with a 
two-dimensional Gaussian, ay = 0.4 fm [18]. The simu¬ 
lation results for central 3 He+Au and d+Au are listed in 
Table [I] with the d +Au values from Ref. [IS’. 


Using track-BBC pairs we construct the distribution 
over relative azimuth, and from that the normalized cor¬ 
relation function: 


S(A<t>,p T ) = 


,/ »rtrack(pT)—PMT\ 

d{ w PMT^ Same event ) 


dA(f> 


C(A^,p T ) = 


S(A(/),pt ) Jo ^ M{Acj>,p T ) dA4> 
M(A</>,p T ) S(Acl),pT)dA(j) 


(1) 

( 2 ) 


The weighting u>pmt for each pair is taken as the PMT 
charge. The signal distribution S is over pairs in the 
same event; the mixed distribution M is over pairs from 
different events in the same event centrality and z ver¬ 
tex bin, and serves to correct for any nonuniformity in 
acceptance over Aip. 



FIG. 1. (Color online) The BBC-S charge distribution in min¬ 
imum bias 3 He+Au events and in high-multiplicity triggered 
events, scaled appropriately. The dashed line indicates the 
threshold selecting the 0%-5% most central events. 


TABLE I. Monte Carlo Glauber characterization results. 


System 

N par t 

IVcoll 

£2 

£3 

0%-5% 3 He+Au 

25.Oil.6 

26.1T2.0 

0.50T0.02 

0.28T0.02 

0%-5% d+Au 

17.3T1.2 

18.1T1.2 

0.54T0.04 

0.19T0.01 


To estimate the contribution to the flow measurements 
from elementary processes, such as jets and resonance de¬ 
cay, we first examine azimuthal correlations in minimum 
bias p+p and central 3 He+Au events across a long range 
in pseudorapidity, with |A?y| ~ 3.5 between tracks in the 
PHENIX central arm at a given pt and charge measured 
in the BBC PMTs. We use the BBCs for these corre¬ 
lation functions because the detector configuration and 
performance is uniform over many years, which enables 
us to combine data for p+p collisions at -/s = 200 GeV 
from running in 2005, 2006, 2008, and 2009. This results 
in 2.7 billion total minimum bias p+p events. 



FIG. 2. (Color online) The azimuthal correlation functions 
C(A(/),pt), as defined in Eq. [5] for track-BBC pairs with 
different track pr selections in (a)-(c) HM 3 He+Au col¬ 
lisions and (d)-(f) minimum bias p+p collisions both at 
y's JVJV = 200 GeV. The track pr bins are (a),(d) 0.4-1.0, 
(b),(e) 1.0-2.0, and (c),(f) 2.0-3.0 GeV/c. Each correlation 
function is fit with a four-term Fourier cosine expansion; the 
individual components n = 1 to n = 4 are drawn on each 
panel, together with the fit function sum. 

Figure [ 2 ] shows the correlation functions C(A<P,pt) 
for different pr bins, for HM 3 He+Au collisions using 
(a)-(c) BBC-S and (d)-(f) minimum bias p+p collisions 
using both BBCs. We analyze these shapes by fitting 
each C(A(/),pt) to a four-term Fourier cosine expansion, 
f{A4>) = 1 + 2 Cu(pt) cos(nA(j)). The sum func¬ 

tion and each individual cosine component are plotted 
in Fig. [ 2 ] for each distribution. Central 3 He+Au col¬ 
lisions show a clearly visible enhancement of near-side 
pairs, producing a local maximum in the distribution at 
A (j) ~ 0. In contrast, p+p collisions can be described 
almost completely by the dipole term cos(A</>), as ex- 
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pected generically from back-to-back jets and transverse 
momentum conservation. 

As in our d+Au analysis | 10 | . we estimate quantita¬ 
tively the correlation strength that would be observed 
in a class of 3 He+Au collisions purely from elementary 
processes, such as in p+p collisions. Approximating the 
3 He+Au collisions as superpositions of some number N 
of p+p collisions, we would then expect the correlation 
strengths c n from the superposition to be the same as 
in p+p, but diluted by a factor of l/N. We then ap¬ 
proximate this dilution factor as simply the ratio of total 
charge observed in the BBC-S detector in p+p versus 
3 He+Au: 


c HeAu elementary _ C P + P {p T ) 


(EQ BBC - S U 

(EQ BBC - S ) H cAu 


(3) 


with a numerical value of 1/(20 .6 ± 0.4). Figure [3](a) 
shows the c n ’s from the 3 He+Au correlation functions, 
and from the p+p with the dilution factor applied. The 
cf 3 is found to be positive and the c^ p to be negative. 
This indicates that the correlation of elementary pro¬ 
cesses contributes positively to the v 2 but negatively to 
the V 3 . The ratios in Fig. |3][b) show that the relative cor¬ 
relation strength in 3 He+Au from elementary processes 
grows with pr, as might be expected from jet processes 
for example, but does not exceed 7% (15%) for c 2 ( 03 ). 

We now quantify the strength of collective behav¬ 
ior through the v n coefficients. The v n coefficients for 
charged hadrons at midrapidity are measured in cen¬ 
tral 3 He+Au events via the event-plane method [25> as 
v n {PT) = (cosn(^> Partlcle (pT) —'k® bs ))/Res('k° bs ), where 
the average is over particles in the same pt bin and events 
of the same centrality. The nth-order event-plane direc¬ 
tion T° bs is determined in each event with the BBC- 
S or FVTX-S detectors. The ’F° bs are corrected for 
each detector with a standard event-plane flattening tech¬ 
nique 1251128] to remove the effect of any small, residual 
nonuniformities in the detector response. As a cross¬ 
check we use event planes from both the full FVTX-S 
covering —3.0 < 77 < —1.0 and a subsection covering 
—2.5 < 77 < —1.5. The choice of the latter is to avoid 
edge effects and still retain good FVTX-S acceptance. 

We calculate the resolutions Res('F(/ bs ) for each de¬ 
tector at each n using the standard three-event-plane 
method uni 1251, combining two event planes with the nth 
order event plane determined from central-arm tracks, 
restricted to low pt (0.2 < pt < 2.0 GeV/c) to mini¬ 
mize contribution from jet fragments. In the case of the 
n = 2 event plane, the resolution is also estimated using 
the first-order event plane measured with spectator neu¬ 
trons in the shower-maximum detector of the zero-degree 
calorimeter 28) [29] on the Au-going side (77 < —6.5). 
The values for the resolution obtained in both methods 
are found to be consistent within uncertainties. 

The event plane resolutions for each detector and order 
are shown in Table [II] The v 2 and V 3 measured using the 
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FIG. 3. (Color online) (a) The c„(pt) coefficients 
for track-BBC pairs from high multiplicity 3 He+Au col¬ 
lisions (filled, denoted “c A ”) and c u {pt) for pairs in 
minimum bias p+p collisions times the dilution factor 
(EQ BBC - S ) p+p /(EQ BBC - S ) HeAu (W d ^oted “c B ” ). 

(b) The ratio c B /c A is shown with statistical error bars. 

(c) Comparison of extracted values of v 2 and V 3 for midrapid¬ 
ity tracks in central 3 He+Au using event plane measurements 
with detectors in different pseudorapidity intervals (see text). 


three event planes described above are shown in Fig.[3](c), 
and are consistent between detectors to within 5%(15%) 
for 7 ^ 2 ( 773 ) over the whole pt range. 


TABLE II. The resolution of nth-order event-plane angles 
measured by the BBC-S and FVTX-S detectors. 


Subsystem 

Res(T? bs ) 

Res(4'? bs ) 

BBC-S (-3.9 < 77 < -3.0) 

0.110 

0.034 

FVTX-S (-2.5 < 77 < -1.5) 

0.232 

0.052 

FVTX-S (-3.0 < 77 < -1.0) 

0.274 

0.070 


The main sources of systematic uncertainty for the v n 
measurements are: ( 1 ) track backgrounds from weak de¬ 
cays and photon conversions; ( 2 ) multiple 3 He+Au col¬ 
lisions in a bunch crossing (pile-up); (3) biases in event 
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plane determination; (4) the effect of detector alignment 
and performance on the v n measurement; and (5) elemen¬ 
tary process/nonflow correlations. We assign the follow¬ 
ing values to account for these systematic uncertainties: 
(1) We estimate the track background contribution by 
reducing the spatial matching windows in PC3 from 3<r 
to 2cr, and find a change of less than 2%(5%) fractionally 
in v 2 (v 3 ). (2) We expect the v n from pile-up events to be 
modestly reduced. Conservatively assuming that pile-up 
events that contaminate the sample at the level of 4%- 
5% have a negligible v n , this results in a systematic 
uncertainty. (3) Event plane effects are estimated from 
v n measurements using different event plane detectors as 
shown in Fig. [3jc) . They are no more than 5%(15%) 
for v 2 (v 3 ). (4) The difference of v n for charged hadrons 
measured by the east and west DC arms are found to 
be less than 2%(15%) for 772 ( 773 ). (5) The contribution 
from nonflow correlations at each pp is estimated from 
Fig- §b), reaching a maximum of 7%(15%) for v 2 (v 3 ). 
We do not attempt to correct for this contribution and 
instead treat it as a systematic uncertainty. All of these 
contributions are summed in quadrature. 

The final v n results are determined using the event 
plane measured in the FVTX-S covering -3.0< 77 <-1.0, 
and these are shown in Fig. [4] with the systematic uncer¬ 
tainties as described above. We observe sizable v 2 and v 3 
anisotropies that both rise as a function of pp- It is no¬ 
table that the v 2 {pt) values for central 3 He+Au collisions 
are very similar within uncertainties with those reported 
earlier in central rf+Au collisions [10]. In scenarios where 
these anisotropies reflect the initial geometry, this sim¬ 
ilarity would be expected as the initial eccentricities e 2 
for central d+Au and 3 He+Au are essentially identical, 
as shown in Table |U The same calculations indicate a 
much larger £3 in 3 He+Au compared with d+Au colli¬ 
sions. However, the d+Au data used in [1(4 were taken 
in 2008, without a central trigger and before the FVTX 
was installed, and did not allow extraction of a statisti¬ 
cally significant v 3 in d+Au. 

We now compare the experimental data with the¬ 
ory predictions in the literature. Four such predictions 
shown in Fig. [4] employ viscous hydrodynamics with 
77 /s at or near the conjectured lower bound 1 / 4-77 30]. 
The Glauber-|-Hydro [3T] (IP-Glasma+Hydro |B2]) utilize 
Glauber (IP-Glasma) initial conditions, and both over 
predict the magnitude of the v 2 and v 3 data. Improved 
agreement may be achieved by utilizing a larger value of 
77 /s or by the inclusion of a transition from QGP to a 
hadronic cascade, which has much larger viscous effects 
and thus decreases the overall flow. The SONIC calcula¬ 
tion m employs Glauber initial conditions, viscous hy¬ 
drodynamics, and then at T = 170 MeV a transition to a 
hadronic cascade. The SUPERSONIC calculation [33] ad¬ 
ditionally includes pre-equilibrium dynamics that boosts 
the initial velocity fields at the earliest times. The im¬ 
pact of pre-equilibrium is modest on the v 2 values and 



FIG. 4. (Color online) Results for v 2 (circles) and 773 (squares) 
as a function of pr for inclusive charged hadrons at midrapid¬ 
ity in 0%-5% central 3 He+Au collisions at y/s NN = 200 GeV; 
error bars are statistical and shaded bars are systematic un¬ 
certainties as described in the text. Also shown are various 
theoretical calculations, see text for details and references. 


both calculations agree with the data within uncertain¬ 
ties. The effect of pre-equilibrium on v 3 is significantly 
larger as the triangular flow takes longer to develop [14]. 
The SUPERSONIC prediction agrees well with the exper¬ 
imental data for p T < 1.5 GeV/c, and then the data 
trends towards the SONIC prediction at higher pp- 

Lastly, we compare to calculations utilizing the A- 
Multi-Phase-Transport (ampt) model [33], which in¬ 
corporates both partonic and hadronic scattering, and 
has recently been compared with anisotropies in cen¬ 
tral p+Pb and d+Au collisions [35] [36]. ampt results 
for 3 He+Au agree reasonably with the experimental v 2 
and v 3 data for pp < 1 GeV and then significantly un¬ 
derpredict the data. Possible underlying causes of the 
anisotropies within the AMPT model are discussed in 
Ref. [37]. 

We have presented first results on azimuthal 
anisotropies v 2 and v 3 in central 3 He+Au collisions at 
y/s NN = 200 GeV. Calculations including hydrodynamic 
expansion of the initial hot spots created in 3 He+Au col¬ 
lisions qualitatively describe the data. Further compar¬ 
ison with different theoretical models should be infor¬ 
mative in terms of the contributions from the initial ge¬ 
ometry and each time stage in the medium evolution in¬ 
cluding pre-equilibrium. Forthcoming results from p+Au 
collisions at RHIC will provide a full suite of geometries 














with highly asymmetric collisions to constrain the origin 
of the observed anisotropies. 
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